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Abstract 

The  discrete  time  detection  of  a known  constant  signal  in  additive  white 
stationary  Laplace  noise  is  considered.  The  receiver  operating  characteristics 
of  the  Neyraan-Pearson  optimal  detector  are  presented  and  compared  with  those  of 
the  linear  detector.  Also,  some  results  obtained  using  a Gaussian  approximation 
to  the  distribution  of  the  test  statistic  are  presented. 


1.  INTRODUCTION 

The  Gaussian  distribution  is  a popular  model  for 
noise  statistics,  and  it  is  valid  in  a wide 
variety  of  situations.  However,  in  some  cases, 
such  as  impulsive  noise,  this  model  may  not  be 
satisfactory.  Impulsive  noise  is  typically 
characterized  as  noise  whose  distribution  has  an 
associated  "heavy  tail"  behavior.  That  is,  the 
probability  density  function  (pdf)  approaches  zero 
more  slowly  than  a Gaussian  pdf.  The  references 
in  [1]  and  (2]  give  a summary  of  some  forms  of 
impulsive  noise  and  situations  where  it  arises. 

Certain  forms  of  impulsive  noise  may  be  charac- 
terized by  a Laplace  distribution.  That  is,  the 
pdf  of  the  noise  is  given  by 

f(n)  = I . (1) 

Notice  that  this  model  has  the  "heavy  tall" 
behavior  associated  with  impulsive  noise.  The 
Laplace  distribution  is  popular  in  statistics  and 
many  of  its  properties  have  been  studied  [3). 
Furthermore,  it  is  used  as  a noise  model  in 
engineering  studies.  For  example.  Miller  and 
Thomas  [1)  used  Laplace  noise  in  a numerical  study 
of  relative  efficiency.  Bernstein,  et.al.  [41 


comment  on  the  non-Gaussian  nature  of  ELF  atmos- 
pheric noise,  and  they  give  a plot  of  a typical 
experimentally  determined  pdf  associated  with  such 
noise  [4,  Figure  10].  This  experimentally  deter- 
mined pdf  Is  similar  to  a Laplace  pdf,  and  on  a 
linear  graph  the  difference  is  barely  distin- 
guishable. Also,  the  limiting  case  of  the  Mertz 
model  (5)  for  the  amplitude  distribution  of  impul- 
sive noise  is  identical  to  the  distribution  of  the 
amplitude  of  Laplace  noise  (6).  Kanefsky  and 
Thomas  [7]  considered  a class  of  generalized 
Gaussian  noises,  obtained  by  generalizing  the 
Gaussian  density  to  obtain  a variable  rate  of  ex- 
ponential decay.  The  J.aplace  distribution  is 
within  this  class  of  generalized  Gaussian  distri- 
butions. Also,  Duttweiler  and  Messerschmitt  [8] 
refer  to  the  Laplace  distribution  as  a model  for 
the  distribution  of  speech. 

In  this  paper  we  are  concerned  with  the  detection 
of  a signal  in  Laplace  noise.  Recently,  we  inves- 
tigated the  distribution  function  of  the  test 
statistic  which  arises  In  this  problem,  and  we 
obtained  a convenient  expression  for  it  [6J.  Here 
we  use  these  results  to  Investigate  various  aspects 
of  the  detection  problem.  In  the  fol lowing  sect  ions 
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wt*  briefly  state  the  problem  a»ul  summarii'e  some 
previous  results.  For  larj*e  sample  sizes*  tlie 
C'entral  l.imii  Thcureri  is  olten  used  to  obtain  an 
ipproximat ion  to  the  distribution  of  the  test 
statistic.  We  compare  the  results  obtained  using 
the  actvjal  distribution  and  using  the  approxi- 
mating Gaussian  distribution.  Then  we  consider 
the  receiver  operating  characteristics  of  the 
Neynian-Pearson  optimal  detector.  The  popular  but 
subopt imai  linear  detector  is  also  considered  and 
is  compared  to  the  optimal  detector. 


2.  I'RELIMINARIKS 

We  consider  testing  for  the  presence  or  absence  of 
a positive,  constant  signal  s in  additive,  sta- 
tionary, white  Laplace  noise.  The  problem  is 
modeled  as  the  following  hypothesis  testing 
problem: 


H 

H 


0* 

1* 


1=1,2 

s > 0 


N 


based  on  N independent  Identically  distributed 
observations,  we  are  to  decide  to  announce  that 
the  signal  is  absent  or  present. 


For  this  problem  both  the  Neyman-Fearson  optimal 
detector  and  the  linear  detector  have  the  fol- 
lowing structure.  The  N observations  are  fed  into 
the  zero  memory  nonlinearity,  g(*),  whereupon  they 
are  summed  to  give  the  test  statlsti<" 

N 

t - [ gCx  ). 

1-1  ' 

The  test  statistic  is  then  compared  to  a threshold 
T from  which  a decision  concerning  the  presence  or 
absence  of  the  signal  Is  made.  The  choice  of  the 
threshold  T determines  the  false  alarm  probability 
(land  the  detection  probability  b.  However,  to 
determine  this  interrelationship,  It  Is  necessary 
to  know  the  probability  distribution  of  the  test 
statistic  t.  This,  In  turn,  depends  on  the  non- 
linearity g(*)»  We  will  now  consider  two  choices 
for  the  nonlinearity  and  will  investigate  the 
corresponding  detector  performance  under  various 
condit  tons. 


3.  THE  NEYMAN-FEARSON  OFTIMAI.  DETECTOR 

By  a Neyman-Fearson  optimal  detector  we  mean  a 
detector  which,  for  a fixed  a,  maximizes  B.  The 
nonlinearity  for  the  optimal  detector  in  our  pro- 
blem is  well  known  19)  and  is  the  amplifier- 
limiter  given  by  the  following  expression: 

I Y S , X > S 

2yx-ys,  0 ^ x < s (2) 

->s,  X < 0 

Consider  first  the  case  where  Hq  is  true;  that  is, 
no  signal  is  present.  If  has  the  Laplace  dis- 
tribution given  in  (1),  then  g(Xj)  can  easily  be 
shown  to  have  the  following  distribution  function: 

(x)=  ^ + expl-  ^ ( tx+v)  J(;(.^y^)dv 


I 


->s 


u(x-Yx)  , 


where  u(x)  denotes  the  unit  step  function  which  is 
one  for  nonnegative  x and  zero  otherwise,  and  G(x) 
denotes  the  gate  function  which  is  one  for 


|.xl  < 1/2  and  zero  otherwise.  Let  F 


(0) 


N 


(•) 


denote  the  distribution  function  of  t for  the 
optimal  detector  under  H . Then 

(x-v)  dFj^°>(v). 

By  a rather  lengthy  but  straightforward  procedure 
(6],  it  can  be  sliown  that* 

.(0) 


N 


(x)  = 2 


[ 

• [x+(N-2p-2q)-ts||j  -u 


exp[-(p+q)>s)'exp[--2  (x+N>s)  | 


(2) 

I x+(N-2p-2q) YS 1 


2 ^ I (m)  ■ u[  x+(N-2m)  TS  ] , 


where  e, ( • ) , 

K 

as 


the  Incomplete 


e^(x) 


k 

I 


m*0 


exponential , 

m 

x 

m! 


is  defined 


We  now  consider  the  signal  present  case,  i.e.  H . 

f 1 1 ^ 

We  let  F,,  ^ (x)  denote  the  distribution  function 

N 

of  the  test  statistic  under  Since  the  Laplace 

pdf  Is  symmetric,  It  can  be  shown  [1]  that 


F<»(x) 


1 


p(0) 

N 


(-X). 


(A) 


*A  typopraphlcnl  error  (omission  of  a pair  of 
brackets)  in  [ft]  has  been  corrected  here. 


Eqs.  (3)  and  (4)  thus  completely  determine  the 
performance  of  the  Neyman-Pearson  opt  Inial  detec  tor . 


i 
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4.  THE  LINEAR  DETECTOR 


By  a linear  detector,  we  mean  a scheme  such  as 
that  described  previouslv,  but  where  t[ie  tunc- 
tion  g(*)  is  g(x)=x.  That  is,  the  test  statistic 
is  simply  the  sum  of  the  observations.  The  linear 
detector  is  Neyman-Pearson  optimal  for  Gaussian 
noise  and  is  a comnionly  used  detector. 


Consider  the  signal  absent  case,  i.e.  H^.  In  this 
situation,  the  test  statistic  is  given  by 


t = 


N 


1 

i=i 


where,  again,  the  are  independent  identically 
distributed  random  variables  with  the  pdf  of 
Equation  (1).  Thus,  the  pdf  of  t,  Ppj(')»  Is 
obtained  from  N-1  convolutions  of  the  Laplace  pdf 
with  itself.  This  is  given  by  [3.  p.  24)  as 


o-(N+k)  (N+k-i)!  (yIx|) 

(N-1)!  k!(N-k-D! 


N-k-1 


After  a straightforward  integration  (10],  we 

obtain  (x) , the  distribution  function  of  the 
N 

test  statistic  of  the  linear  detector  under 


N-1 


k»0 


,-(N+k) 


/N+k-l\ 

( ^ 


c4°>(x) 


■‘’N-k-l^^*^  , X i 0 


<r(-x). 


x < 0 


(5) 


In  the  signal  present  case,  the  tesk  statistic  is 

given  by  N 

t - l X,  + Ns  , 

1-1  ‘ 

where,  once  again,  the  are  Independent  iden- 
tically distributed  random  variables  with  the 
density  function  of  Eq.(l).  Let  G^^^(x)  denote 
the  distribution  function  of  the  test  statistic 
of  the  linear  detector  under  H^.  Then  we  have 


c‘*’(x)  - C.|J°Nx  - Ns)  . 


(6) 


Equations  (5)  and  (6)  completely  determine  the 
performance  cf  the  linear  detector. 


5.  THE  GAUSSIAN  APPROX IM/‘ i iON 


In  non-Gaus**lan  detection  problems  of  the  type 
considered  in  this  paper,  the  derivation  of  an 


expression  for  the  distribution  function  of  the 
test  statistic  for  the  Neyman-Pearson  optimal 
detector  is  frequently  a mathematically  intractable 
problem.  In  many  such  cases,  for  sufficiently 
large  N,  an  appeal  is  made  to  the  Central  Limit 
Theorem  to  arrive  at  an  approximation  for  the  dis- 
tribution function  of  the  test  statistic.  Thus  it 
is  instructive  in  the  present  case  to  compare  the 
exact  results  with  those  resulting  from  the 
Gaussian  approximation. 


Let  X be  a random  variable  with  the  density  func- 
tion of  Eq.(l).  Let  g(*)  be  the  optimal  non- 
linearity given  by  Eq.(2).  Then 

E{g(X))  - / g(x)  1 dx  . 

A straightforward  integration  yields 


E(g(X)  ] = 1 - ys  - e ■ 

Similarly,  we  get 

VAR|g(X)l  - /lg(x)  - 1 + YS  + e e dx 

= 3 - 2e-^=  - Ays  . 

Thus  the  mean  and  variance  of  t under  Hq  are, 
respect ively , 

E^lt)  ” Nil  - ys  - e = m 

VARjjlt]  = N^[3  - 2e"''®  - Aysc"^®  - . 

Using  the  relation  in  Eq.(4),  It  follows  that  the 
corresponding  values  under  are  given  by 


Ej^l  t } “ ^ 

VARjIt]  -=  VAR^lt]  * . 

Let  Ij^^\x)  and  denote,  respectively,  the 

Gaussian  approximations  to  the  distribution  func- 
tions of  the  lest  statistic  under  and  Hp  Then 


and 


i*'bx)  - , 

N 0 


where 


i(x) 


-v^/2 


dv  . 


Let  and  respectively,  denote  the  false 
alarm  probability  and  the  detection  prob«ability 
resulting  t rom  the  Gaussian  assumption.  Then  we 
have 


I 

i 

,t 

I 

j 


! 


2 


(7) 


U 


made  using  the  results  presented  in  the  earlier 
sections. 


and 


. 1 - >(M) 


, . . 


In  praciue»  one  may  use  Eq.(7)  to  set  the  value 
of  the  threshold  F.  For  example,  if  ‘»s=l,  N=15, 
and  the  desired  false  alarm  probability  is  0.3, 
the  Gaussian  approximation  yields  a threshold  of 
approximately  1.208  and  a detection  probability  of 
approximately  0.b2H.  however,  using  Eqs.(3)  and 
(4)  we  find  that  for  this  threshold,  a 0.02 
ind  P 0.41.  In  fact,  for  u=  0.3,  we  find  that 
the  actual  detection  probability  is  greater  than 

0.49.  Thus,  In  this  case,  the  Gaussian  approxi- 
mation is  extremely  conservative.  In  Table  I we 
•.  .'mpare  the  actual  values  of  i and  8 for  the 
optimal  detector  against  and  8^  for  several 
values  of  I when  N*25  and  rs»0.5.  It  is  seen 
t rom  the  table  that,  in  this  case,  the  Gaussian 
approximation  is  not  verv  good  (even  though  N»25). 


1 

1 

B 

4.032 

0.002 

0.278 

0.278 

0.452 

i.024 

0.007 

0.308 

0.443 

0.487 

1.512 

0.035 

0.356 

0.696 

0.540 

0.000 

0. 123 

0.407 

0.877 

0.593 

-1.512 

0.  104 

0.460 

0.965 

0.644 

-2.016 

0.384 

0.477 

0.978 

0.660 

Table  1.  Exact  values  of  a and  B ^nd  those 

resulting  from  the  Gaussian  approxi- 
mation, for  several  values  of  the 
threshold;  ys=0.5,  N=25. 
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Figure  1.  Performance  of  Optimal  Detector  for 
>s»0.3  and  Different  Values  of  N. 


Figure  2.  Perfornvmce  of  Optimal  Detector  for 
N“10  and  Different  Values  ot  >s. 
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Figure  3.  Performance  of  Optimal  Detector 

Ci'mpared  to  Performance  t'f  Linear 
Detector,  N-H),  s“l,  v“0.5. 


Figure  4.  Performance  of  Optimal  Detector 

Ci'inpared  to  Performance  of  Linear 
Detector,  N“20,  s»0.1,  >■!. 


